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PREFACE 


Since  the  promulgation  of  the  General  Theory  of  Relativity, 
science  has  come  to  look  upon  gravitational  phenomena  as  a  direct 
consequence  of  the  curvature  of  a  four— dimensional  space— time 
continuum  arising  from  a  momentum-energy  distribution.  While 
this  viewpoint  enables  one  to  predict  observed  results  with  greater 
accuracy  and  is  more  satisfying  in  certain  basic  aspects  than  the 
Newtonian  formulation,  it  requires  solutions  of  a  comparatively 
complicated  system  of  nonlinear  partial  differential  equations  of 
the  second  order.  In  an  attempt  to  shed  some  light  on  the  properties 
of  the  solution  manifold  of  the  general  relativity  field  equations  with 
arbitrary  momentum-energy  tensors,  we  presented  analyses  of 
certain  aspects  of  the  general  second— order  discontinuity  problem. 
The  underlying  idea  behind  this  approach  is  that  the  analysis  of 
discontinuity  properties  are  usually  considerably  simpler  than  the 
analysis  of  the  complete  field  equations — particularly  if  the  field 
equations  are  nonlinear.  Following  the  ideas  delineated  in  our 
previous  work  concerning  the  evolutionary  properties  of  the  second- 
order  problem,  we  present  in  this  Memorandum  a  general  dynamical 
theory  of  discontinuity  surfaces  and  the  associated  jump  strengths  of 
both  physical  and  geometrical  quantities.  The  results  reported  here 
form  the  basis  for  a  general  analysis  of  galactic  morphology  which 
will  be  presented  in  succeeding  communications. 
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SUMMARY 


Let  E  be  a  time-like  hypersurface  in  an  Einstein— Riemann 
space  E  which  is  a  support  hypersurface  for  jumps  in  the  mo¬ 
mentum-energy  tensor  and  in  at  least  one  second  coordinate  deriva¬ 
tive  of  the  metric  tensor  h^g(x).  A  system  of  necessary  conditions 
is  obtained  for  the  existence  of  solutions  of  the  Einstein  field  equations 
in  the  presence  of  jumps  on  E.  These  conditions  are  shown  to  be 
expressible  in  terms  of  surface  tensors  and  surface  tensorial  dif¬ 
ferential  systems.  In  particular,  a  system  of  algebraic  surface 
tensor  equations  are  obtained  for  the  jump  strengths  of  the  metric 
tensor  and  a  system  of  first  order  covariant  surface  differential 
equations  are  obtained  for  the  jump  Strengths  of  the  momentum- 
energy  tensor.  These  equations  involve  only  the  surface  components 
of  the  jump  strengths  and  the  components  of  the  first  fundamental 
form.  An  additional  algebraic  condition  is  obtained  which  involves 
the  surface  jump  strengths  of  the  momentum— energy  tensor  and  the 
second  fundamental  form.  If  it  is  assumed  that  one  knows  the  mo¬ 
mentum-energy  tensor  on  one  side  of  E,  explicit  formulae  are 
obtained  for  the  continuation  of  the  momentum-energy  tensor  across 
the  discontinuity  hypersurface  E  by  integrating  the  differential 
system  governing  by  the  jump  strengths  of  the  momentum— energy 
tensor.  This  integration  is  effected  by  means  of  two  principles 
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which  state  a  fundamental  relation  between  geometry  and  physics 

and  the  fact  that  physical  systems  are  governed  by  second  order 

differential  relations.  The  case  in  which  the  hypersurface  E  forms 

# 

a  static  three-dimensional  hyperbolic— normal  metric  space  E  is 
then  examined.  Results  are  thus  obtained  which  will  be  of  signifi¬ 
cance  in  the  application  of  the  discontinuity  method  to  problems  in 
relativistic  cosmology. 
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1.  INTRODUCTION 

Since  the  promulgation  of  the  General  Theory  of  Relativity, 
science  has  come  to  look  upon  gravitational  phenomena  as  a  direct 
consequence  of  the  curvature  of  a  four-dimensional  space-^time 
continuum  arising  from  a  momentum— energy  distribution.  While 
this  viewpoint  enables  one  to  predict  observed  results  with  greater 
accuracy  and  is  more  satisfying  in  certain  basic  aspects  than  the 
Newtonian  formulation,  it  requires  solutions  of  a  comparatively 
complicated  system  of  simultaneous  nonlinear  partial  differential 
equations  of  the  second  order.  Some  success  has  been  achieved  in 
obtaining  solutions,  but  the  methods  employed  rely  heavily  on  a 
large  number  of  assumptions,  and  usually  involve  very  restricted 
forms  of  the  momentum-energy  tensor  or  ignore  this  tensor  alto¬ 
gether.  With  the  exception  of  existence  theorems  for  the  initial- 
value  problem,  very  few  results  have  been  obtained  with  respect  to 
_ 

The  RAND  Corporation,  Santa  Monica,  California 
** 

University  of  California — San  Diego,  La  Jolla,  California 


-2- 


the  properties  of  the  solution  manifold  with  arbitrary  momentum- 
energy  tensors. 

In  an  attempt  to  shed  some  light  on  the  properties  of  the  so¬ 
lution  manifold  of  the  general  relativity  field  equations  with  arbitrary 
momentum-energy  tensors,  we  presented  analyses  of  certain  aspects 
of  the  general  second-order  discontinuity  problem  [  1  through  6]  . 

The  underlying  idea  behind  this  approach  is  that  the  analysis  of  dis¬ 
continuity  properties  are  usually  considerably  simpler  than  the 
analysis  of  the  complete  field  equations — particularly  if  the  field 
equations  are  nonlinear.  Such  analyses  also  provide  the  basis  upon 
which  significant  information  concerning  the  detailed  physical  pro¬ 
cess  can  be  gleaned;  see,  for  example  [7],  [8]  and  [9] .  Following  the 
ideas  delineated  in  our  previous  papers  concerning  the  evolutionary 
properties  of  the  second-order  problem,  we  have  obtained  a  general 
dynamical  theory  of  discontinuity  surfaces  and  the  associated  jump 
strengths  of  both  physical  and  geometrical  quantities.  These  results 
are  presented  here  and  will  form  the  basis  for  a  general  analysis  of 
galactic  structures. 

2.  PRELIMINARY  CONSIDERATIONS 

In  this  section  we  collect  certain  required  results  from  the 
differential  geometry  of  hypersurfaces.  We  state  these  without  proof 
and  refer  the  reader  to  the  standard  texts  or  to  the  exposition  given 
in  [  6  ]  . 
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Consider  a  metric  space  of  the  type  used  in  general  relativity 
theory;  that  is,  an  Einstein-Riemann  space  E  whose  metric 
structure  is  defined  by  the  quadratic  differential  form 

(2.1)  ds2  =  hAB  dxA  dxB,  (A,  B  =  0,  1,  2,  3), 

having  signature  —2  and  coefficients  h^B  which  are  functions  of 
A 

the  coordinates  x  of  the  space  E.  As  indicated  in  (2. 1),  capital 
Latin  indices  will  have  the  range  0,  1,  2,  3  and  will  be  summed 
over  this  range  in  accordance  with  the  summation  convention.  Let 
E  denote  a  regular  hyper  surface  in  E;  that  is,  E  can  be  defined 
parametrically  by  equations 

in  n\  A  ,A.  0  12. 

(2.  2)  x  «  f  (u  ,  u  ,  u  ), 

where  the  fA(u°)  are  continuously  differentiable  functions  of  the 

Of 

parametric  or  surface  coordinates  u  such  that  the  functional 
A  ct 

matrix  ((8f  / 3u  ))  has  rank  3  for  all  values  of  the  u's  under 
consideration.  Here  and  throughout  this  paper  lower  case  Greek 
indices  will  be  associated  with  surface  quantities  and  obey  the 
summation  convention  with  the  range  0,  1,  2. 

Over  such  a  hypersurface  E  one  can  define  a  normal  vector 
with  covariant  components  by  the  equations 


_v  A  . .  n  .  A  def 
(2.3)  *,Na  =  0  <*,  > 


8^(10/811“). 
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We  shall  assume  throughout  this  paper  that  E  is  a  time— like 
hypersurface,  and  hence  its  normal  vector  is  space  like.  The 
normal  vector  to  E  may  thus  be  normalized  by  the  requirement 

(2.4)  NaNA  =  -1. 

A 

The  quantities  x^  for  a  =  0,  1,  2  are  the  components  of  three 
independent  contravariant  vectors  in  E  at  any  point  of  E  and  can 
also  be  interpreted,  for  A  fixed,  as  the  components  of  four  co- 
variant  vectors  on  the  surface  E. 

The  coefficients  a^  of  the  first  fundamental  form  of  the 

(X 

surface  E  are  given,  as  functions  of  the  surface  coordinates  u  , 
by  the  equations 


(2-5)  ■  hAB  \  */3 


As  defined  by  (2.  5),  the  quantities  &q^  transform  according  to  the 

tensor  law  when  the  surface  coordinates  undergo  their  admissible 

transformations.  The  surface  tensor  determined  by  the  a^  is 

called  the  metric  tensor  of  the  surface  E.  It  can  be  shown  that  the 

first  fundamental  form  of  E  is  nonsingular  and  that  it  has  signature 

aB 

— 1.  Also,  denoting  the  inverse  of  a^  by  a  ,  it  can  be  shown 
that 


c.  aB  A  B  ,  AB  A  .TB 
(2.6)  a  x^  X£  *  h  +N  N  . 
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Before  proceeding  further,  it  appears  advisable  to  state  the 
exact  assumptions  of  continuity  and  differentiability  which  will  be 
involved  in  this  paper.  For  this  purpose  let  us  denote  by  J  the 
region  (open  set)  consisting  of  some  neighborhood  of  E  which 
contains  £  in  its  interior.  Denote  by  the  subregion  of  J  lying 
on  one  side  of  £  and  by  J ^  the  subregion  lying  on  the  other  side. 
Let  us,  furthermore,  denote  by  the  domain  +  £  (that  is, 
the  point  set  consisting  of  the  subregion  and  the  points  of  the 
surface  £  as  boundary  points),  the  domain  J,  +  £  being  denoted 
by  D„.  The  following  assumptions  are  now  made. 

u 

A  ot 

Aj!  The  functions  f  (u  )  in  equations  (2.2)  defining  the 

3 

hypersurface  £  are  of  class  C  . 

A  :  The  metric  components  h^g(x  )  are  functions  of  class 

1  3 

C  in  the  region  J  and  of  class  C  in  the  domains  and  D£- 


It  follows  from  these  assumptions  that  the  Christoffel  symbols 
A 

of  the  Einstein— Riemann  space  E  are  continuous  across  £. 
Also,  from  assumption  A^  and  equations  (2.6)  it  is  seen  that  the 
metric  components  a^  (u  )  are  continuous  on  £  and  have  con¬ 
tinuous  first  partial  derivatives;  hence  the  Christoffel  symbols 

CK 

Apy  determined  by  the  quantities  a^  are  continuous  functions  of 
the  surface  coordinates.  One  can  thus  construct  the  first  surface 


covariant  derivatives  of  differentiable  tensorial  quantities  defined  on 
E.  Thus,  in  particular,  we  have 
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N 


a;l 3 
A 


=  x 


•a 


N 


Q,P 

A 


A  aV 
—  x  A'a 
y  ap 


+  x 


B 


.A 


aA  C 

+  N  abc  Xo 


for  the  components  of  the  covariant  derivatives  of  the  mixed  surface 

A. 

and  space  vectors  defined  on  E  by  x^  and  the  unit  normal  vector 
N  respectively.  Use  of  the  semicolon  to  denote  covariant  differ¬ 
entiation  will  be  continued  throughout  this  paper;  correspondingly  we 

shall  use  the  comma  to  denote  partial  differentiation  with  respect  to 

A  a 

the  space  coordinates  x  or  the  surface  coordinates  u  as  in¬ 
dicated  by  the  indices. 

Denoting  the  coefficients  of  the  second  fundamental  form  on 
E  by  b^  we  have 

,  A  XT 

b  „  =  — x  a  N  .  . 
afi  or;p  A 


Hence  the  b  „  are  continuous  functions  of  the  surface 

aP 

and  are  moreover  symmetric  in  the  indices  a  and  P 

A 

quantities  xa.p  are  symmetric  in  these  indices.  The 
b occur  in  the  following  important  relations 


(2-7>  Vs  *  b«s  «A’ 


(2.8)  N 


.  Py  A 
=  b  „  or'  X  . 


;a  aP 


coordinates 
since  the 
functions 
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3.  JUMP  CONDITIONS 

So  far  the  hypersurfaoe  £  has  been  any  regular,  time— like 
hypersurface  in  a  four-dimensional  Einstein-rRiemann  space  E.  We 
now  restrict  our  attention  to  those  hypersurfaces  in  E  which  carry 
basic  field-theoretic  information  in  the  sense  that  they  are  the  support 
hypersurfaces  for  the  field  discontinuities.  In  view  of  our  previous 
assumptions,  this  is  accomplished  by  the  following  requirement. 

A0:  There  is  a  discontinuity  in  at  least  one  of  the  second 

u 

derivatives  of  the  functions  h^^(x  )  at  points  of  the  hypersurface  £. 


i  A.  .  .  I 

The  symbol  |W  *  '  *  |  will  be  used  to  denote  the  jumps  in  the 

D  •  •  • 

A 

quantities  W  *  '  ’  across  Z;  that  is.  the  differences  in  the  limits 

Di  •  « 

A. .  . 

of  W  *  *  '  as  points  in  Z  are  approached  from  J-  and  J0 
respectively.  By  assumption  A^,  we  have 


(3.1)  lhABl  =0,  IHABiCl  =0, 


while  Ag  states  that  CP  I  ^  ®  for  some  choice  of  the  indices. 

In  fact,  it  can  be  shown  [  6]  that  there  exist  functions  X.^gfu**)  of 
the  coordinates  of  the  hypersurface  E  such  that 


(3.2) 


AB,CD' 


XAB  NC  ND‘ 


In  view  of  (3.1)  and  the  fact  that 


the  X.'s  may  be  viewed  as  the  jump  strengths  of  the  h— field. 


We  now  impose  the  requirements  that  the  structure  of  the  space 
E  is  determined  by  the  Einstein  field  equations 


<3-4>  BAB"IBhAB 


K  T 


AB‘ 


Here  B^^  are  the  components  of  the  Ricci  tensor  of  E,  B  is  the 
scalar  curvature  B 


AB 


AB 

h  ,  are  the  components  of  the  momentum- 


energy  tensor,  and  k  is  the  so  called  gravitational  constant.  Let 
us  denote  the  jump  strengths  of  the  momentum— energy  tensor  across 
L  by  S^.  that  is 


(3.5) 


SAB  1 TAB 1 » 


ot 

where  the  S's  are  functions  of  the  surface  coordinates  u  .  It  is 
then  easily  shown  (see  [  2  ] )  that  a  necessary  condition  for  the 
existence  of  solutions  to  the  Einstein  field  equations  under  assumptions 
A  ,  A  ,  A  is  that  the  discontinuity  strengths  satisfy  the  equations 

(3-6>  ^AB1  'I  |E|  hAB  '  *  SAB' 

The  system  (3.6)  also  may  be  viewed  as  conditions  for  the  con¬ 
tinuation  of  solutions  to  the  Einstein  equations  across  surfaces  of 
discontinuity  of  the  momentum-energy  tensor  (see  [4  ] ). 


Denoting  the  components  of  the  complete  curvature  tensor  of 
A 

E  by  B  BCD,  we  have 


_.A  .A  .A 

R  =  A  —  A  +  # 

BCD  BC,D  BD,  C  ’ 


where  the  asterist  denotes  terms  which  are  quadratic  in  the 
Christoffel  symbols.  When  use  is  made  of  these  relations  and 
equations  (3.1)  and  (3.2)  we  obtain 


i  A  i  Am 

(3.7)  2|BBCD|.h  kMCNBND+XBDNMNc-VBCNMND 


-xmdnbV> 


from  which  the  left-hand  side  of  (3.6)  may  be  evaluated.  Introducing 
the  quantities 

(3.8)  <j>A  =  XabNB,  «|>  =*aNA, 


the  explicit  evaluation  of  (3.6)  is  found  to  be  given  by 


<3'9)  VAB  +  2  *<A  V  -  +  NA  NB 

CD 

-«.  +  XCDh  )hAB*2,,SAB- 

g 

If  we  multiply  both  sides  of  (3.  9)  by  N  ,  sum  on  the  repeated 
index  B,  and  use  equations  (3.8),  the  left  members  of  the  resulting 
equations  vanish  identically.  We  are  thus  left  with  the  simple  set  of 


relations 
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(3.10)  SabNB  =  0. 


It  is  also  evident  that  the  functions  <f>A  are  undetermined  by  the 
equations  (3.8),  (3.9)  and  hence  may  be  taken  as  arbitrary  functions 

Qf 

of  the  coordinates  u  on  E.  These  functions  clearly  reflect  the 
arbitrainess  in  the  choice  of  coordinates  of  the  space  E .  In  addition, 
if  there  are  no  physical  jumps  (i.e.,  =  0),  the  admissible  choice 

<}>A  =  0  anihilates  the  \'s.  In  this  sense,  the  metrical  jump  strengths 
determined  by  the  <f>'s  have  no  intrinsic  physical  meaning. 

The  fundamental  system  of  equations  (3.  9)  and  (3. 10)  involve 
the  jump  strengths  \  .  _  and  S._,  which  are  tensor  quantities  under 
admissible  coordinate  transformations  of  the  four-dimensional  space 
E.  Now  consider  the  quantities  as  defined  by  the  equations 

<3-“»  s«h‘sab\  ■ 

(X 

Under  admissible  transformations  of  the  surface  coordinates  u  , 
these  quantities  transform  according  to  the  indicated  tensor  law, 
and  hence  constitute  the  components  of  a  surface  tensor.  A  direct 

calculation  based  on  (2.6)  and  (3. 10)  gives 

. .  _AB  _or/3  A  B 

(3. 12)  S  *  S  p  xa  , 

aB 

where  the  S  are  obtained  from  when  we  raise  the  indices 

aB 

by  means  of  the  contravariant  quantities  a  in  the  usual  manner . 
Equivalently,  if  we  assume  (3.12),  equations  (3.10)  are  identically 
satisfied.  We  thus  see  that  the  surface  quantities  SQ^  give  a  unique 
determination  of  the  components  over  any  specified  surface  E, 

and  that  this  determination  identically  satisfies  equations  (3.10). 

One  may  also  verify  the  following  result: 
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(3. 13)  S  ■  SAB  hAB  -  Sab<*A  *=  3*  -  NA  NB) 
‘  <SAB  \  ^  '  S<,0  a"S  • 


We  now  consider  the  surface  tensor  whose  components  are 
defined  by  the  equations 


(3. 14)  X  a  =  X  xA  x? 

ap  AB  a  p 


Let  us  write  (3.9)  in  the  equivalent  form 


(3. 15)  1AB  -  2  «  Sab  -  2  *(A  Nb)  -  <♦  ♦  «  S)  Na  Nb 


-  *  S  hAB 

A  B 

and  contract  both  sides  with  x^  x^  .  The  result  is  the  following 
simple  system  of  equations: 


(3.16)  ltf  -S.jS). 


It  is  easily  seen,  when  use  is  made  of  (3. 12)  and  (3.16),  that  the 
relations  (3. 15)  can  be  written  in  the  equivalent  form 


(3.17)  XAB  =  2(X°^  -  X  a“^)  xA  xB  -  2  4>(A  NB) 

X)  NA  NB  +  X  hAB, 


where  X  is  defined  by 
(3.18)  X  *  \Qp  . 
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qB 

Conversely  the  equations  (3.17)  in  which  the  X  are  derived  from 
(3. 16)  imply  the  relations  (3.  9)  on  E. 

Combining  the  considerations  of  the  previous  paragraphs,  we 
are  thus  led  to  the  following  basic  result.  Necessary  conditions  for 
the  existence  of  solutions  to  the  Einstein  field  equations  under  as¬ 
sumptions  A^,  A^,  Ag  are  given  by  the  tensorial  surface  equations 


(3.19)  X 


*(2S«e  "SV- 


in  which  the  surface  quantities  and  X.^  are  determined  from 

the  jump  strengths  SAB  and  XAB  by  (3.11)  and  (3. 14)  respectively. 
By  this  result  we  are  permitted  to  base  our  succeeding  considerations 
on  surface  tensors  and  the  intrinsic  geometry  of  E. 


4.  DIFFERENTIAL  RELATIONS 

It  is  seen  from  the  results  of  Sec.  3  that  knowledge  of  the 

surface  quantities  is  sufficient  to  determine  the  X^  and 

hence  the  X  to  within  the  geometrical  structure  of  E,  i.e.  to 
AB 

A 

within  the  functions  a^  ,  NA>  and  hAB  on  E.  As  yet,  however, 
we  have  not  used  the  full  content  of  the  Einstein  theory  for  the 
relations 

(4.1)  T®.b  •  0 

are  still  at  our  disposal.  We  shall  now  show  that  the  system  (4. 1) 
leads  to  a  set  of  differential  relations  on  E  which  will  partially 
determine  the  structure  of  and  the  geometry  of  E. 


We  first  note  that  from  (3.5)  we  have 


(4.2)  s®  =  It® I  =  It®  I  xc . 

A ;a  A  ;a  A;C  a 

If  we  multiply  (4.  2)  by  x®  a°^,  sum  on  the  repeated  index  a  and 
use  (2.  6),  we  are  led  to  the  equations 


<4'3»  SA;«  ^lT^cl<hCD+NCND>. 


Obvious  manipulations  of  (4.3)  then  give 


(4.4) 


ap  D,  „B.t 
“  ^hCD-FANC- 


where 

<4'5»  FA"|TA;C|NC 

are  functions  of  the  coordinates  u*  of  the  hypersurface  E.  Now, 
since  (4. 1)  holds  on  both  sides  of  E,  the  tensor  equations 


<4-6>  ^Ib1'0 

must  hold  on  E.  Substituting  from  (4.4)  into  (4.6),  we  then  obtain 
the  differential  relations 


(4.7) 


SAB;a 


B 

X0 


where  the  quantities  F^,  which  are  given  by 
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constitute  the  components  of  a  spacial  vector  on  the  hypersurface  L. 

We  now  replace  the  differential  relations  (4.  7)  by  a  set  of 

equations  which  involve  the  quantities  rather  than  the  in 

accordance  with  the  viewpoint  expressed  at  the  end  of  Sec.  3.  For 

A 

this  purpose  let  us  multiply  both  sides  of  (4.  7)  by  N  to  obtain 


(4.9) 


_tA  0  B  oP  _tA  def 

N  S .  _  x0  a  =  N  F  .  =  v 

AB;a  0  A  * 


Then  by  covariant  surface  differentiation  of  (3. 10)  and  use  of  (2. 8) 
we  find  that 


(4. 10)  b-  +  X  =  0. 


On  the  other  hand,  if  we  contract  both  sides  of  (4.  7)  with  x^,  we 
have 

(4. 11)  S  x®  xA  a°^  =  F  xA  d=f  F  . 

AB  ;a  p  y  A  y  y 

It  now  follows  from  (2.  7),  (3. 10)  and  (4.  11)  that 


„a  .BA  afi.  B  A  a0  afi.  B  A. 

Sr,o  '  ®AB  Xe  \  a  >„  '  SAB;«  \  *  +  SAB  *  (x0  Xy>;« 

•  F  +  S  a°^(NB  b  .  xA  +  xf  NA  b  )  •  F  . 
y  AB  ap  y  p  ay  y 


We  have  thus  proved  the  following  fundamental  result.  Necessary 
conditions  for  the  existence  of  solutions  to  the  Einste  n  field  equations 
under  assumptions  A^,  A^,  A^  are  given  by  the  relations 
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(4.12)  S“S  bafS  +  x  -  0, 

(4.13) 

over  the  surface  £,  where 

(4. 14)  X  -  Fa  NA,  F„  ■  Fa  xA  ,  Fa  -  ITa.c  I  NC  Nb. 

Although  we  will  not  pursue  this  approach,  the  above  results 
may  be  used  to  obtain  equations  for  the  surface  quantities  X.^. 

To  see  this,  we  first  solve  (3. 19)  for  s<*  and  thereby  obtain 


(4.15)  2  k  =  \ap  \  aa£  . 


Substitution  of  (4. 15)  into  (4. 12)  and  (4. 13)  then  leads  to  the  desired 
results,  namely 


(4.16)  \a<3  b^  -2  U1  +2«x  =  o, 


(4.17) 


\/3  =  2kF(S> 


where  ft  =  (1/2)  b  0  a  is  the  mean  curvature  of  L  in  E. 

ap 


5.  THE  CONTINUATION  PROBLEM 

The  results  established  in  the  previous  section  show  that  the 
jump  strengths  of  the  momentum— energy  tensor  and  its  covariant 
derivative  cannot  be  specified  in  an  arbitrary  fashion  in  view  of  the 
conditions  imposed  by  the  Einstein  field  equations.  However,  a 
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certain  degree  of  freedom  is  allowed.  To  see  this,  we  first  note 
that  the  F  depends  on  the  S's  and  the  geometry  of  £  in  accordance 
with  the  equations 


(5.1) 


_A  _a 3  ,  A  ,  _a£  A 
F  =  S  b  o  N  +  S  a  x 
ap  ;P  a 


This  follows  from  (4. 12),  (4. 13),  (4. 14)  and  the  fact  that  (N^)  and 

(x  ;  for  a  =0,  1,2  form  a  vector  basis  on  E.  It  is  thus  evident 
a  ’  ’ 

that  if  we  specify  S<*P  in  an  arbitrary  fashion  and  use  our  previous 
results,  the  quantities  \q^  and  F^  are  uniquely  determined  in  a 
manner  consistent  with  the  Einstein  field  equations. 

Suppose,  now,  that  we  know  the  momentum-energy  tensor  on 
one  side  of  £  and  hence  the  limiting  values  of  this  tensor  and  its 
divergence  as  E  is  approached  from  that  side.  The  question  then 
arises  as  to  the  continuation  of  the  momentum-energy  tensor  across 
E  such  that  the  Einstein  field  equations  will  be  soluable.  Since  we 
have  implicitly  assumed  that  we  do  not  know  T^  on  the  other  side 
of  £,  we  cannot  use  the  field  equations  to  determine  the  h's  and 
hence  the  X's.  We  can  therefore  not  use  (3. 16)  to  determine  the 
S's  since  we  do  not  know  the  X,'s.  It  thus  follows  that  the  only 
information  at  our  disposal  is  provided  by  the  existence  conditions 


(5.2) 


and 
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(5.3)  b^  +  y  =0. 

As  previously  noted,  we  could  specify  the  S's  and  then 

A 

determine  the  quantities  F  by  (5.  1).  This  procedure  would 

amount,  however,  to  fitting  the  equations  of  motion  (i.e.,  the 

divergence  of  the  momentum-energy  tensor)  on  one  side  of  E 

to  those  on  the  other  side  by  a  specification  of  the  quantities 

|T®  |  N  ,  these  latter  quantities  being  uniquely  determined 

AjC  d 

by  the  arbitrarily  choses  S^.  Conceptually,  equations  of  motion 

are  in  some  respects  more  fundamental  than  momentum— energy; 

B  C 

in  addition  the  quantities  |T.  |  N  N  determine  certain 

A  jL/  o 

properties  of  the  motion  normal  to  E  and  hence  are  more  easily 
observed  than  the  S's.  Hence,  having  observed  that  such  an 
artificial  fitting  process  is  possible,  we  shall  disregard  it  as  a 
general  method  of  procedure.  We  are  thus  left  with  the  problem 
of  determining  the  S's  such  that  equations  (5.  2)  and  (5.  3)  are 
satisfied  when  the  functions  and  ^  are  arbitrarily  pre— 

Ct 

assigned  functions  of  the  coordinates  u  . 

Let  us  first  fix  our  attention  on  the  system  (5.  2).  These 
equations  may  be  considered  as  an  invariant  differential  system 
in  the  three-dimensional,  hyperbolic-normal  metric  space  E* 

a 

with  coordinates  u  and  metric  differential  form 

(5.4)  da ^  =  a^(u^)  du*  du^  . 
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Because  of  the  three-dimensional  character  of  E*,  the  components 
K a@y\  curvature  tensor  can  be  expressed  by 


(K«-IKa«), 


where  are  the  components  of  the  permutation  tensor  of  weight 

X  aB 

zero,  Kap  =  ^  Qp\>  anc*  K  =  a  .  The  Einstein  tensor  on 

E*  thus  serves  to  determine  the  complete  curvature  tensor 

If  the  system  (5.2)  is  to  possess  a  solution  in  E*,  there  must 

aB 


exist  a  tensor  Q 


such  that 


^aB  ,„a  _  /3a . 

(5.5)  Q  =  F  ,  (F  =  a  ). 

aB 

In  the  context  of  differential  equations,  the  tensor  Q  plays  the 
role  of  a  particular  solution  of  (5.  2).  The  fact  that  (5.2)  is  a  linear 

aB 

differential  system  in  S  allows  us  to  write  the  general  solution 
in  the  form 


(5.6)  =  Q0f!  +  Zafi  , 


0(3 

where  Z  is  the  general  solution  of  the  homogeneous  system 


(5.7)  =0. 


To  some  extent  it  would  appear  that  we  have  robbed  Peter  to 
pay  Paul,  for  we  now  have  to  face  the  problem  of  determining  the 
functions  .  Now,  the  system  (5.  7)  is  formally  similar  to  the 
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equations  of  motion.  In  addition,  the  T's  represent  the  momentum- 
energy  complex  in  E  while  the  S's  represent  part  of  the  jump 
strength  of  momentum-energy  and  hence  have  energetic  interpreta¬ 
tions  in  E*.  The  formal  analogy  goes  even  deeper,  however,  for 
the  Z's  also  partially  determine  the  \'s  and  these  latter  quantities 
represent  the  jumps  in  the  second  coordinate  derivatives  of  the 
geometrical  quantities  h^g.  Also,  as  is  well  known,  the  intrinsic 
geometry  of  a  discontinuity  hypersurface  enters  into  the  determi¬ 
nation  of  classical  jump  strengths  in  an  essential  fashion.  It  is  thus 

natural  to  assume  that  the  intrinsic  geometry  of  E*  partially 

aB 

determine  the  S's.  Now,  the  quantities  Q  are  a  particular 

solution  of  (5.  2)  and  hence  depend  on  the  specific  physical  quantities 

Fa.  Equations  (5.  6)  thus  show  that  the  dependence  of  Sa ^  on  the 

intrinsic  geometry  of  E*  (that  is,  on  functions  of  K  ^  ^  )  must 

aB 

arise  through  the  quantities  Z  .We  are  therefore  led  to  state  the 
following  postulates  for  the  integration  of  the  system  (5.7). 

aB 

P^:  The  quantities  Z  are  the  components  of  a  metric  tensor 
differential  invariant  of  the  space  E*. 

P2:  This  metric  tensor  differential  invariant  is  of  the  second 
order  and  linear  in  the  second  derivatives  of  the  components  of  the 
metric  tensor  of  E*. 

Postulate  P2  expresses  the  usual  restriction  to  second  order 
relations  which  is  commonly  assumed  in  physical  theories.  On  the 
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other  hand,  postulate  P1  together  with  (5.6)  is  a  fundamental 
statement  concerning  the  relation  between  physics  and  geometry. 

The  procedure  is  now  straightforward.  It  follows  from  (5.  7) 

that  the  divergence  of  the  metric  tensor  differential  invariant  (Z°^) 
must  vanish  identically.  From  the  known  procedures  for  con¬ 
structing  such  tensor  differential  invariants  [10],  it  follows  that  the 
most  general  quantity  which  satisfies  our  requirements  is  given  by 

(5.9)  Zep  -  8  (KaS  -  f(K  +  n)ao(J), 

where  ©  and  n  are  constants.  "We  shall  henceforth  assume  that 

these  constants  have  fixed  values,  and  in  particular  that  ©  ^  0. 

We  have  now  established  a  consistent  procedure  whereby 

S's  may  be  obtained  such  that  we  simultaneously  satisfy  the  existence 
requirements  (5.2),  for  (Fq)  an  arbitrarily  assigned  surface  vector, 
and  the  requirements  of  postulates  P^  and  We  still  have  to 
satisfy  the  existence  requirement  (5.  3).  Now,  the  Q's  area 

particular  solution  of  the  linear  system  (5.  2)  and  hence  we  may  write 

(5.10)  ^  +  Pa(3  , 

where  the  W's  are  unique  functionals  of  the  F's  such  that 

(5.11)  ■  Waf(XFy> 

is  an  identity  in  X.  and  the  P's  are  any  functions  of  the  coordinates 
u*  for  which 
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(5.12)  P®^  0  -  0. 

;p 

If  we  substitute  (5.6)  and  (5.10)  into  (5.3),  we  obtain 

(5.13)  Pa\p  =  -  X  “(Wa/3  +  Za\p  . 

Since  all  of  the  terms  on  the  right-hand  side  of  this  equation  are 
known,  (5.13)  is  seen  to  be  an  equation  for  the  determination  of 
Combining  this  with  (5. 12),  we  have  four  equations  for  the 
determination  of  the  six  quantities  Hence  the  existence 

requirements  (5.  3)  can  always  be  satisfied;  in  fact,  we  are  free  to 
add  two  more  conditions  in  most  cases. 

Combining  the  above  results,  we  have  the  following  con— 

or 

elusion.  Let  the  functions  ^  &ud  F  be  arbitrarily  assigned 

a 

functions  of  the  coordinates  u  .  Then  a  consistent  procedure  for 
the  continuation  of  T^g  across  E,  such  that  (1)  Postulates  Pj 
and  P2  are  satisfied  and  (2)  the  Einstein  field  equations  are  soluble, 
is  given  by 

(5‘14)  S«|3  +  ®(K«0  “|<K+  n>V’ 

(5.15)  Ka/3b^  ~  0(k+  n)n  +  Qa\p  +  y  -  o, 

(5.!6)  Kafi  =2«<0Ka(3  +Qa0)-|(©(K-  n)  +  2Q)a^, 


(5.17) 
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where  ®  and  n  are  constants. 

6.  THE  DISCONTINUITY  HYPERSURFACE  AS  A  STATIC  3-SURFACE 

The  problem  of  continuing  the  momentum— energy  tensor  across 
L  arose  because  of  an  assumed  lack  of  knowledge  of  T^g  on  one 
side  of  L.  There  are  a  number  of  important  physical  problems  in 
which  this  is  the  case  and,  in  addition,  the  functions  h^.g(x  )  are 
unknown — either  because  the  Einstein  field  equations  have  not  been 
solved  or  because  there  is  insufficient  information  to  effect  such  a 
solution.  Fortunately,  this  lack  of  information  is  usually  compen¬ 
sated  by  known  properties  of  the  discontinuity  hypersurface.  We 
shall  now  consider  a  particularly  important  example  of  this  situation. 

Let  L  be  a  hypersurface  for  which  the  three-dimensional 
metric  space  £*  is  known  to  be  static;  that  is,  the  space  E*  admits 
a  one-parameter  group  of  isometries  whose  trajectories  form  a 
time— like  normal  congruence  [11].  This  means  that  there  is  a  vector 
field  Y  defined  on  L*  such  that 


(6.1) 


Y(«;0)  =  °»  Y 


a  2^ 

r  Yq  ,  ■  0,  Y  Y  «  e  . 

l«  rJ  '  <*  1 


Q 

where  ^  is  a  finite-valued  function  of  the  coordinates  u  .  From 
the  physical  point  of  view,  these  conditions  state  that  the  infinitesi¬ 
mal  2-spaces  of  observers  traveling  along  the  curves  of  the  con¬ 
gruence  (i.e.,  the  curves  defined  by  du°/dp  =  YQ(vP)  where  p  is 
the  parameter  on  the  congruence)  mesh  into  finite,  space— like 
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two-surfaces  called  space  sections.  These  space  sections  are 
mapped  isometrically  onto  each  other  by  the  group  of  point  trans¬ 
formations  generated  by  the  congruence.  The  space  sections  of 
E*,  and  hence  of  E,  correspond  to  what  we  would  normally  refer 
to  as  the  two-dimensional  boundary  of  a  three-dimensional  space- 

like  body  obtained  from  E  by  an  appropriate  time— section.  In 

.A.  ct  A 

addition,  the  vector  field  with  components  Y  =  Y  in  E  is 

everywhere  tangent  to  E  and  transports  the  first  fundamental  form 

on  E  without  change.  Hence,  the  requirement  that  E*  be  static 

is  equivalent  to  the  requirement  that  E  be  generated  by  the  two- 

dimensional  boundary  of  a  body  acted  upon  by  a  one-parameter 

group  of  point  transformations  which  preserves  the  metric  structure 

on  that  is  is  metrically  stable  in  E . 

It  should  be  noted  that  the  requirement  that  E*  be  static  does 

not  necessarily  imply  that  the  Einstein-Riemann  space  E  is  static. 

For  instance,  suppose  that  there  is  a  C*  vector  field  defined 

—  A 

in  E  and  such  that  X.x  =Y  ,  where  the  bar  denotes  evaluation 

A  a  a’ 

on  E.  We  then  have 


Y.  =  (X.  xA)  a  =X.  /QxA.  +X.  a 
(«;£)  A  {a) ;0  A;(/3  a)  A  ap 

-X/A  d.xAx?  +  X  a  a  . 

(A;B)  a  p  A  aP 


Hence,  if  E*  is  static,  can  vanish  for  bQ^  £  0  only  if 

—  A 

X^N  vanishes. 
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In  order  to  obtain  the  implications  of  the  system  (6. 1),  we 
define  a  vector  U  in  L*  by  the  equations 


(6.  2)  U  =  e  ®  Y 
a  a 


Using  (6.  2}  to  eliminate  the  vector  Y  from  the  system  (6. 1),  we 
obtain  the  equivalent  system  (see  [12]  ,  pp.  65  f) 


(6.3)  U  U"  =  1,  U  a  =  fr  U„,  =  U  a  =  ¥  . 

a  ’  a;p  a  p’  a  a;p  .a 


As  a  direct  consequence  of  this  system  we  have 


(6.4)  2  U  '  r  n  .  U0*2ft  ,  Ufll  +  2  fj  A  a  U  ■,  . 

<*;[0y]  «;[y  Pi  «  [0  y] 

Hence,  if  we  contract  (6.4)  on  the  indices  (a,  y)  and  note  that  (6.  3) 
implies 


ti  U«  =  -  tf  ft  Ufl  , 
a;P  o  P’ 

we  finally  obtain  the  relations, 

(6.5)  K  „  U*  =  (a"y*  )  Ufl  . 

op  ;ay  p 

The  content  of  equations  (6.5)  is  that  the  Ricci  tensor  of  L*  admits 

U  as  a  time— like  eigenvector  and  that  the  corresponding  eigenvalue 

aB 

is  the  D'Alembertan  of  i.e.,  a  ¥  Q ^  . 

We  now  apply  these  conditions  to  the  continuation  problem  for 
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T  .  If  we  contract  (5.14),  we  have 
A  B 

(6.6)  ©  K  =  -2(S  -Q  +  |  CD  n)  , 


where  we  have  used  the  obvious  notation 
Qap  a* ^  .  We  may  then  solve  (5.14)  for 


S  =  Sap  a.a ^  and  Q  = 
K  „  to  obtain 


(6.7)  »Ka/1  -S^  -(S-Q  +  . 


Using  (6.  7)  and  the  fact  that  has  to  satisfy  the 

conditions  stated  by  (6.5)  if  £*  is  static,  we  can  obtain  the  corre¬ 
sponding  requirements  which  must  satisfy.  These  requirements 

are 

(6.8)  (Sag  -Q„s)  ■  «>  Ua  , 

for  p  given  by 

(6.9)  p  *  ©  a"^  +  S -Q  +®n  . 

Equation  (6. 8)  states  that  admits  U  as  a  time— like 

eigenvector  with  associated  eigenvalue  p;  assuming  that  this  eigen¬ 
value  is  simple,  we  obtain  a  representation  for  S^,  namely 


*«„(!  +0U„  u0  +  V 

The  tensor  with  components  which  appears  in  (6. 10)  is 

arbitrary  to  within  the  requirements  that  it  be  symmetric  and  admit 


-26- 


U  as  a  null  vector.  The  representation  (6. 10)  gives 

(6.11)  S=Q+p+cr“, 

and  hence  (6.9)  leads  to  the  interesting  result. 

(6.12)  0  (a®^  +  n)+a®  -  0. 

Since  we  have  exact  requirements  of  the  U's,  namely 
equations  (6.3),  we  may  obtain  significantly  more  information  con¬ 
cerning  the  structure  of  for  static  spaces  L*.  Without  loss 

of  generality,  we  may  assume  that 


(6.13)  -  p  U„  Ufl  +Q„ft  +M 


a 


where  the  quantities  p  and  are  to  be  determined  so  that  U 

is  an  eigenvector  of  M.  Noting  that  (6.3)  implies 


(6.14)  *  U®  »  0,  if  =  0, 

f 

the  substitution  of  (6. 13)  into  the  equations  (5.  2)  leads  to 

(6. 15)  i  U  +  n  tT  +  a  =  0, 
r  a  r  a  Qf;p  ’ 

and  hence 

ji  +  a  if  =0. 

^  a;P 


If  we  use  the  fact  that 


and  then  add  and  subtract  the 
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quantity  (n^)  a  ,  the  system  (6.15)  becomes 

(6.16)  (Ufl  +(Mf  +li®fif)^  =0. 

Thus,  since 

(6.17)  V  Dy)(6j  -Up  U\ 

and 

(6.18)  IVT  U  =  4  UH 

we  have  the  seven  equations  (6. 12),  (6.17)  and  (6. 18)  for  the 
determination  of  the  eight  functions  4,  (A,  and  . 

As  an  indication  of  the  use  of  these  results,  consider  the 
particular  case  in  which  the  function  p  is  constant.  The  system 
(6. 16)  then  reduces  to 

(6.19)  (M^  +  =0, 

and  hence  it  is  possible  to  consider  the  particular  situation  in  which 
we  have 

(6.  20)  . 

The  equations  (6.13)  and  (6.17)  then  give 

<6'21>  S«g  +  >*<Ua  U0  “*V> 


and 
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(6.22)  -Uff  U^. 

Hence,  substituting  (6.  22)  into  (6. 17)  we  obtain  the  important  result 

(6.23)  0(a“^  +  n)  —  2  p.  =  0. 

The  function  V  thus  satisfies  a  Schroedinger  type  equation,  and 
in  addition,  must  be  such  that 

(6.24)  Ua  =  0. 

This  last  requirement  follows  from  (6.  14). 

Although  the  above  example  is  highly  artificial  in  the  manner 
in  which  it  has  been  introduced  here,  related  situations  arise  in 
certain  fundamental  problems  associated  with  relativistic  cosmology 
and  galactic  structure.  These  problems  will  be  treated  in  suc¬ 
ceeding  papers  in  which  equations  similar  to  (6.  23)  lead  to  new  and 


fundamental  results. 


-29- 


REFERENCES 


1.  Thomas,  T.  Y. ,  "On  the  Propagation  and  Decay  of  Gravitational 

Waves,"  J.  Math.  Anal,  and  Applications,  3(1961),  pp.  315— 

335. 

2.  Edelen,  D.  G.  B.,  and  T.  Y.  Thomas,  "Discontinuities  in  the 

Einstein  Field  for  General  Momentum— Energy  Tensors,  "  Arch. 
Rational  Mech.  Anal. ,  9  (1962),  pp.  153—171. 

3.  - ,  "Differential  Compatibility  Conditions  on  the 

Momentum-Energy  Tensor  and  Necessary  Conditions  for  the 
Existence  of  Solutions  to  the  Einstein  Field  Equations, "  Arch. 
Rational  Mech.  Anal. ,  9  (1962),  pp.  245—254. 

4.  Edelen,  D.  G.  B.,  "On  the  Continuation  of  Orthogonal  Structure 

Across  Surfaces  of  Discontinuity  in  the  Momentum— Energy 
Tensor,"  J.  Math.  Anal,  and  Applications,  4  (1962),  pp.  346— 
351. 

5.  Thomas,  T.  Y. ,  "Gravitation  and  the  Origin  of  Material 

Particles, "  Arch.  Rational  Mech.  Anal.,  10  (1962),  pp.  197— 

204. 

6.  - ,  "Hypersurfaces  in  Einstein-Riemann  Space  and 

their  Compatibility  Conditions,  "  J.  Math.  Anal,  and  Applications, 
(this  issue). 

7.  Coburn,  N.,  "Discontinuity  Relations  for  Charged,  Com¬ 

pressible,  Relativistic  Fluids,"  J.  Math,  and  Mechanics, 

10  (1961),  pp.  361-391. 

8.  Israel,  W.,  "Discontinuities  in  Spherically  Symmetric 

Gravitational  Fields  and  Shells  of  Radiation,  "  Royal  Soc. 

(London),  Proceedings,  Series  A,  248  (1958),  pp.  404-414. 

9.  Synge,  J.  L.,  Relativity:  The  General  Theory,  North— Holland 

Publishing  Co. ,  Amsterdam  (1960). 

10.  Thomas,  T.  Y. ,  Differential  Invariants  of  Generalized  Spaces, 

Cambridge  University  Press,  London  (1934). 

11.  Lichnerowicz,  A.,  Theories  Relativistes  de  la  Gravitation  et 

de  L'Electromagnetisme,  Masson  et  Cie,  Paris  (1955). 

12.  Witten,  L.,  Gravitation:  An  Introduction  to  Current  Research 

(more  laterj!  "  ——  ~~ 


